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Spatial behavior of spin transport in a Sierpinski gasket fractal is studied from two dimensions 
to quasi-one dimension subject to the Rashba spin-orbital coupling. With two normal metal leads 
represented by self-energy matrix, discretizing the derived continuous Hamiltonian to a tight-binding 
version, Landauer-Keldysh formalism for nonequilibrium transport can be applied. It was observed 
that the spin Hall effect presents in the distribution of spin density critically depends on the fractal 
structure and the shape of the thinfilm. The local spin density and transmission are numerically 
tested by the present quantum transport calculation for the fractal and non-fractal thinfilm varying 
from two-dimensional square- lattice into quasi-one dimensional fractal shape (Sierpinski triangles). 
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I. INTRODUCTION 

Intensive efforts on spin Hall effect (SHE) both experi- 
mentally and theoretically have successfully built another 
milestone in condensed matter physics^. Spin separa- 
tion in semiconductors is not only possible but natural, 
so that manipulating spin properties of charge carriers 
in electronics is promising 2 ^—. Recently different types 
of SHE were reported on different geometrical setup of 
both magnetic^ and non magnetic materials^. Local 
spin Hall effect often attributed as extrinsic SHE^~— , and 
intrinsic SHEJ^. The intrinsic SHE, spin separation is 
solely due to the underlying spin-orbit coupling in the 
band structure, so that SHE can exist even in systems 
free of scattering within a finite size syste m 13 i 14 . How- 
ever, experimentally, most observations so far have been 
attributed to the extrinsic SHE. More recently nonlocal 
spin Hall effect was also found in various systems ^ 11 d 5 , 
including graphene ^ 16 ' 17 . Considerable effort in this di- 
rection has already revealed the unique features of spin- 
polarized electron transport in the so-called two-terminal 
nano- or mesoscopic devices. Indeed, this is a fast devel- 
oping field, and has also stimulated lot of theoretical work 
based on the non-equilibrium Green's function approach 
within the density functional theory-i£~— . 

A two-terminal device is essentially a single path de- 
vice, SHE studies on a two-terminal device with defects 
show a lot interesting results and the local spin den- 
sities can be affected seriously even by a small point 
defect a 20 i 21 for the breaking of translational symmetry. 
In the present studies we undertake an in-depth study of 
the spin transport in a fractal network based on the non- 
equilibrium Green's function approach, the Landauer- 
Keldysh formalism (LKF) 1 - 7 . - — . We choose a fractal ge- 
ometry following a Sierpinski gasket (SPG) 22 i 23 and ex- 
amples of planar gasket is shown in Fig. [TJ With the 
self-similarity of the system and the enlargement of the 
fractal structure, our studies provide an understanding 
of the spin transport behavior transition from two di- 
mensional system to quasi-one dimensional system. Our 
motivation for the investigation of spin transport within 
fractal system is mainly resulted from our previous stud- 



ies on impurity system. The defect and impurities seri- 
ous affect the local charge density and local spin density 
for the influences on spin transport, moreover, the rela- 
tions between spin current and charge current within a 
film with defect behave very different from that within 
an ideal fil m 20 i 21 . Even though the point-defect induced 
translational symmetry breaking and interference among 
different conducting channels provide subtle reasons for 
the variation of the spin density and current density, the 
detailed analysis of defect and impurities of the multi- 
terminal systems is very rare and the multiply connected 
fractal geometry is not even been probed. The state-of- 
art lithographic technologies made all artificial structures 
become possible. Therefore, the studies of spin transport 
of self-similarity of systems such as the SPG provide the 
possibility of investigating the local density of charge and 
spin on different multiply connected fractal geometry. 

Secondly, a very unique phenonemon in spin transport 
remains to be understood is the spin processional phase 
through different travelling channel, therefore, the quan- 
tum interference among the confined paths of the fractal 
structures will be an interesting and provide a good way 
of studying the self-similarity with reducing scale. Even 
though the charge transport on percolating clusters had 
been studied but the local charge density and local spin 
density, and the Aharonov-Bohm effect in a fractal sys- 
tem without translational invariance remains unclear till 
no w 24 ' 25 . Moreover, the density of transistors on an in- 
tegrated circuit doubles every couple of years. The scale 
of devices is one of important factors. Fractal structures 
also can study the limits of the self-similarity with reduc- 
ing scale for a profound 2D-system. 

The paper is organized as follows, we present 
Landauer-Keldysh formalism with a tight-binding frame- 
work to model the fractal system in Section II. Numerical 
results of local charge density and local spin density of 
different size of fractal system arc studied and compar- 
ison among continuous film in low dimensional system 
also been discussed in Section III. We draw our conclu- 
sions in Section V. 
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II. NON-EQUILIBRIUM GREEN FUNCTION IN 
SIERPINSKI GASKET 



A. The Sierpinski gasket 

The Sierpinski gasket is named after the Polish mathe- 
matician Sierpinski in 1915^ and is one of the well-known 
fractal geometry. Since fractal geometry was established 
in 1967, there are several ways to obtain Sierpinski gas- 
ket (SPG)SiSi22r^Si2L2£. Following by previous studies, we 
choose SPG as our sample and four generations of SPG 
are shown in Fig. [TJ The SPG is perfectly self similar, 
an attribute of many fractal structures. Any triangular 
portion is an exact replica of the whole structure (Fig. 
[[]). The dimension of the SPG is log3/log2 = 1.5849 and 
lies dimensionally between a line and a plane. 
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(a) First generation of SPG. (b) Second generation of SPG. 
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(c) Third generation of SPG. (d) Fourth generation of SPG. 

FIG. 1: (Color online) Fractal structures of SPGs: an 
example of planar gasket with the self-similarity of the 
system. Four SPGs are attached to two semi-infinite 
metallic leads, left lead and right lead, on the edges (yellow 
lines). All SPGS are 32 x 32 square lattice and the white 
regions are the atomic sites. The widths connects with the 
each two leads are (a)16a in first generation of SPG, (b)8a 
in second generation of SPG, (c)4a in third generation of 
SPG and (d)2a in fourth generation of SPG. 



B. The Landauer Keldysh formalism 

The Landauer Keldysh formalism is also called non- 
equilibrium Green's function (NEGF) formalism^ and 
in principle carries all the physical information of the 
investigated under a non-equilibrium but approximately 
steady process^. It is usually used to simulate current 



and charge density in atomic-scale quantum mechanical 
transport when bias applied. Also it can be used to cal- 
culate spin density and transmission in different genera- 
tions of sample. We describe the conductor by an array 
of numbered lattice sites, n = 1,2,..., A x N. For a 
32 x 32 sample contains a square lattice of 1024 sites, as 
shown in Fig. [TJ The corresponding Hamiltonian matrix 
is written in the site basis and each element of matrix 
H m<n records the relation between sites m and n. After 
we label the sites of lattice, we will show that the in- 
teraction between the sites by using the nearest neighbor 
hopping to describe the interaction, viz, the tight-binding 
model. With the help of the second quantization, the 
one-particle Hamiltonian composes of on-site terms and 
a hopping term, 



£ n, C 



(i) 



where cjj(c„) is the electron creation (annihilation) op- 
erator, e n is the on-site energy at n-th site, and t m .„ 
is the hopping matrix. Using the system of the Rashba 
spin-orbital interaction (RSI)2I for example, the hopping 
matrix is t m>n = —tol — it so (a x e n _>. m ) • e z where I is the 
2x2 identity matrix, to is the kinetic hopping strength 
and t so is the hopping strength of RSI. Applying the near- 
est hopping only, one writes 



H 



m.n 



c} l (-toI—tt S o((TXe n ^, ni )-e z 



(2) 



where (m, n) means m — th site interacts with the near- 
est sites n. Further, |r m — r n | = a with lattice spac- 
ing a is satisfied. Since the Hamiltonain has established, 
then considering the ideal (free of spin-orbital orbit in- 
teractions) leads are connecting with the sample left and 
right, respectively. Deriving the self-energy from RefJ^, 
we obtain the self-energy matrix of the left(right) lead, 
^left(right)- Next, using both the tight-binding Hamilto- 
nian and lead self-matrices, one can build up the retarded 
Green function matrix 



G R (E) = [EI — H — Z left {E) - V right {E)]- 



(3) 



where the I is 2A x 2A identity matrix, H is the 
tight-binding Hamiltonian matrix of the conductor, and 
^ieft(right)(E) is the self-energy matrix of the left(right) 
lead. Next, the lesser self-energy matrix is given by 

£<(£) = - ]T [Z p (E-eV p )-i:l(E-eV p )]f (E-eV p ) 

p=L,R 

(4) 

where /o is the Fermi function and eV p is the electric bias 
potential energy applied on lead p = {L(left), R(right)}. 
In the numerical calculation, we set eV/ e /t = +(eVo)/2 
and eV r ight = — (eVo)/2. Moreover ,we consider the zero 
temperature simplistically in order to ignore thermal fluc- 
tuation. Thus the Fermi function reduces to step func- 
tion. Finally, we can get the lesser Green function via 
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the kinetic equation 

G<(E) = G R {E)^<{E)G A {E) (5) 

where G A = {G R } f Thus we get the lesser Green's func- 
tion G < which can be used to calculate the local charge 
density, 

/oo 
dETr s [G< n ] (6) 
-oo 

where the trace calculates in spin space and N n is the 
electron number operator at site n. Also the local spin 
density will be 

fe 1 />oo 

(5 " )= 2 2™ J_ dETr ^°nJ (7) 

here S n = (S x , S y ,S z ) is spin operator at site n of the 
conductor. We assume our devices lay on x and y direc- 
tions. Thus, the current goes along in x-direction. The 
study of the latter, the local spin density, can be cat- 
egorized into the out-plane ((S z ) — (Sj_)) and in-plane 
((S x ,S y ) = (Sy)) components. And we set the hopping 
energy to = 1, the mass of electron m e = 1 , the charge of 
electron 5 = 1, and the Plank's constant divided by 27r, 
= 1, as the units. We also can calculate transmission, 

T = Tr s [T left G R T right G A ] (8) 

where T left{right) = t[T, left(right) - ^ left(right) ] describes 
the coupling of the conductor to the two leads. 

III. NUMERICAL RESULTS AND 
DISCUSSIONS 

A. First Generation of SPG 



When raising into high bias, electrons with high energy 
injected from the lead and transport through the device. 
Now high energy electrons can go from the linear region 
into the QD-like region and then back to the linear re- 
gion again (Fig. [3ajl . Therefore, charge did not localized 
in the QD-like region as in Fig. [5a] and conducting path 
penetrating into the QD-like region, i.e., the bottom of 
the device. Comparing with the low bias situation, both 
the down-spin and up-spin electrons can form conducting 
channels. While within the QD-like region, there are still 
mainly up-spin electrons and only a few down-spin elec- 
trons trapped in the center of the bottom of the device 
at high bias (Fig. [3b|) 

In first generation of SPG, quantum dot-like region oc- 
curs at low bias. Beside the bottom square shape of 
the first generation of SPG is QD-like and spin-polarized 
electrons driven by the RSI will let the up-spin electrons 
trapped within the QD-like regions. Therefore, a spin- 
dependent potential also form within the QD region and 
thus enhance the confinemnt of up-spin electrons. There- 
fore, first generation of SPG behaves as good spin filter 
with RSI. 
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(a) The local charge density (b) S z 

FIG. 2: (Color online) Local spin density and charge 
density under a low bias in first generation. In Fig. I2bl 
there is a linear wire interacting with a quantum dot. 



The first generation of SPG is setting in the band bot- 
tom Eb — — 4io and the Fermi- level Ef = — 3.8ir> Here, 
the black regions of the sample are empty sites and the 
white regions are atomic sites for transporting charges 
and spins, as shown in Fig. [5J The width connects with 
the each two leads is 16a. 

From Fig. [2 under low bias situation, the charge trans- 
port separated into two channels and the local charge 
accumulate in the bottom channel. The first generation 
of SPG look like a combination of linear conducting line 
coupled with a localized quantum dot(QD). A QD is a 
part of material confined in three dimension c 32 ' 33 . Fig. 
I2bl clearly indicates the up-spin electrons mainly local- 
ized at QD-like region and the down-spin electrons is the 
conducting carrier for first generation of SPG. The sep- 
aration of up-spin electrons and down-spin electrons are 
resulted from the coupling of RSI. Therefore, first gener- 
ation of SPG behave as good spin filter under RSI. 
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(a) The local charge density (b) S z 

FIG. 3: (Color online) Local spin density and charge 
density under a strong bias in first generation. 
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B. Second Generation of SPG 

The second generation of SPG (Fig. Ilbl) is composed of 
three self-similar building blocks of first generation (Fig. 
Ta| . There are empty sites (black regions) and atomic 
sites (white regions) in the device. The conditions are 
setting in the band bottom Eb = — 4to and the Fermi- 
level Ef = — 3.8io, as the same as first generation. How- 
ever, the connecting width with two leads is different and 
is 8a instead of 16a. 

The second generation of SPG geometrically look like 
three first generation of SPG connects in series (Fig. [Ta| 
and Fig. Ilb[) . However for low bias situation, the elec- 
trons are mainly accumulated in the bottom of the device 
in second generation of SPG (Fig. 0a|. There is almost 
no clear conducting channel can be observed with the low 
bias. Additionally, comparing with the first generation of 
SPG, the second of SPG is topological different from the 
first generation of SPG with a hole at the center of fractal 
structure. As reported in ReL^S, a point-defect induced 
symmetry breaking significantly changes the conduction 
channel. The charge density seems to be in favor to pass 
through the bottom of the device, instead of through the 
upper line. The width of the top region of the device 
can change the conducting path of devic o 20 i 21 . The lo- 
cal charge density indicated the charge accumulates in 
corners of the building blocks in Fig. [3a] Local spin 
density also indicated that the spin actually associated 
with charge and accumulated at the same sites but with 
alternating polarization of the spin (Fig. I4bl) . Never- 
theless, both local charge density and local spin density 
show that the second generation of SPG is much local- 
ized than the first generation of SPG. A possible reason is 
from the interference of spin-polarized electrons through 
the channels^ around the central holes and thus the con- 
ducting paths almost disappear in low bias. It should be 
noted that there is also no clearly QD-like behavior in 
second generation. Focusing on S z , the 2 nd generation 
of SPG can be interpreted as the superposition of three 
building blocks of I s * generation of SPG (Fig. [2b] and 
Fig. I2a[) but with one bottom down-spin block and two 
up-spin blocks atop (Fig. libl) . 

As raising higher bias, QD-like behavior appear at the 
bottom building block of the device (Fig. I5b[) but similar 
QD-like behavior does not occur in the top two build- 
ing blocks. The condition of confinement is highly re- 
quired without electron leaking. Here a higher bias did 
not provide enough energy for the electrons hopping into 
the bottom building block and also jumping out again. 
Therefore, a QD-like behavior only observed in bottom 
building block. Even though the geometrical similarity 
within SPG, the electrodes only connects to the left and 
right ends of the upper building blocks and thus the sim- 
ilarity of three building blocks was removed by the bias. 
Therefore, the top two blocks become conducting chan- 
nels and the bottom one acts as QD-like region. We also 
noted that the Fermi wavelength and the unit length of 
building blocks are also critical for the conducting chan- 



nels, As the conducting distance between two leads equals 
to the multiple of Fermi wavelength, several oscillating 
periods of local spin density and local charge density can 
be form in SPG^. 
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(a) The local charge density 
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FIG. 4: (Color online) Local spin density and charge 
density under a low bias in second generation. 
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(a) The local charge density 
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FIG. 5: (Color online) Local spin density and charge 
density under a strong bias in second generation. 



C. Third Generation of SPG 

The 3 rd generation is setting in the band bottom 
E b = — 4t and the Fermi-level Ef = -3.39t . Here 
the sites of the lattice are reduced much more (Fig. flc)) 
and the width connecting with two leads become 4a, so 
that the Fermi-level increases. In 3 rd generation, the 
central defect enlarges and there are also other smaller 
defects. Therefore there are many vortexes around the 
building blocks. Spin and charge are accumulated at the 
corners of building blocks. Even though the geometrical 
self-similarity of the 3 rd generation of SPG, due to the 
complicated paths, the transports of self-similarity of the 
I s * and 2 nd generation of SPG are hardly to see even un- 
der high bias in Fig. I7bl Under low bias, the charge 
density appears only in the left of the device and it in- 
dicates that the bias is weak enough to drive the charge 
to the other end of SPG for too many empty sites (Fig. 
aj. For a high bias, local charge density is around the 
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empty sites, however, it does not directly go along the 
narrow line to the other end (Fig. [Taj. It also noted 
that the local spin density shows only the up-spin elec- 
trons were selected to transport within this downward 
triangular SPG (Fig. [7b]). 
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(a) The local charge density 
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FIG. 6: (Color online) Local spin density and charge 
density under a low bias and Ef = —3.69 in third generation. 
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(a) The local charge density 
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FIG. 7: (Color online) Local spin density and charge 
density under a strong bias and Ef = —3.39 in third 
generation. 



D. Fourth Generation of SPG 



The fourth generation is setting in the band bottom 
Ei, = —Ato- Because reducing the sites of the sample and 
the width, 2a, connecting with two leads, the Fermi-level 
moves to — 2.9io- Under the low bias, the spins and the 
charges transport hardly in the region of the upper line 
in Fig. [5J The self-similar feature of the fractal is not 
observed at all. Under high bias, there are accumulated 
charges and spins in the corners of building blocks of the 
sample in Fig. |Hl Similar selection of spin orientation 
was observed in this downward triangular SPG. 
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(a) Local charge density 
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FIG. 8: (Color online) Local spin density and charge 
density under a low bias and Ef — —2.9 in fourth generation. 
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FIG. 9: (Color online) Local spin density and charge 
density under a strong bias and Ef = —2.9 in fourth 
generation. 



E. Transmission 

In a 2D square sample, the transmission starts from 
— 4io shown as blue line in Fig. HOal For 1 st generation 
of SPG (brown) and 2 nd generation of SPG (purple), the 
transmissions are almost the same as in square sample 
(Fig. HOap but with curve slope change to linear depen- 
dence of energy. The lower threshold of transmission of 
3 rd generation of SPG shifts from — 4io to — 3.39to in Fig. 
II Obi and the transmission of 4 th generation shifts further 
to —2.9to- Our results indicated that the transmissions 
of fractal structure are affected both by the shape of the 
sample and the defects within the samples seriously. For 
1 st generation of SPG, the leads contact the sample with 
16 sites (16a) on both left and right. For 2 nd generation, 
the left and right leads contact the sample only with 8 
sites (8a). Since the width of two leads are compara- 
tively wider than the other higher generation of SPG, 
and thus the spin-polarized electrons can be through two 
leads smoothly. For 3 rd generation of SPG, the two leads 
are in contact with only 4 sites (4a) in each side and the 
transmission reduces heavily. It should be noted that the 
central defects of the higher generation of SPG make their 
transport behavior is topological different from the lower 
generation of SPG. The threshold of the transmission 
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FIG. 10: (Color online) Transmission of SPG versus bias 
energy. The bias energy is normalized with the hopping 
energy to. (a) The transmission of a 2D square reference 
sample (blue line) starts from around — 4to and then 
increases with energy. The non-zero transmission of higher 
generation of SPG shifts to smaller values, (b) Zoom in on 
the higher generation of SPGs with the central defects and 
the transmission was suppressed significantly. Note that the 
threshold of higher generation of SPGs also shifts to left and 
the central triangle defects are with effect in third generation 
(green dagger-line) and fourth generation (red triangle-line). 



shifts from — 4io to — 2Mq with increasing number of de- 
fects, i.e., a higher generation of SPG. For the 4 th gener- 
ation of SPG, a visible gap even appears within — 2.7io to 
— 2.5to (Fig. UObp and this arises from the self-similarity 
of SPGs. For the higher generation of SPG, the contact 
sites with leads reduces and the transmission also reduces 
accordingly for the narrowing of the conducting channel 



and most charges and spins are trapped at the bottom of 
fractal structure. A self-similarity behavior of transmis- 
sion can be observed from 2 nd generation of SPG for the 
topological difference with the lower generation of SPG 
(Fig. ESQ. 

IV. CONCLUSION 

In conclusion, we had studied spin transport in a pla- 
nar Sierpinski gasket fractal formed by elementary trian- 
gle unit. We also compare the local charge density, the 
local spin density and also spin Hall effect on the pla- 
nar Sierpinski gasket with a continuous planar and lin- 
ear structure. For four different planar Sierpinski gasket 
fractal structures in our analysis, with a low bias, there 
arc little fractal-like patterns in first generation and sec- 
ond generation. Moreover, because of too many holes on 
the third generation and fourth generation, the electrons 
become very localized and transport can be only observed 
along the linear edge. Therefore, hardly any fractal pat- 
terns were observed for both the local spin density and 
local charge density. Indeed there is quantum dot region 
was observed in first generation (Fig. [2] and Fig. [3]) and 
also barely was observed in second generation (Fig. [5]). 
The behavior can be considered to be a similar arrange- 
ment of single level QD a 36 i 37 sitting at the vertices of a 
continuous line. Also the spin density and the charge den- 
sity indicate that the transmissions in fractal structures 
are smaller than linear structures. And this is consistent 
with the results of transmission. A very interesting phe- 
nomenon was observed in our studies, due to the fractal 
shape, spin transport can be separated from the charge 
transport and details mechanism still need to be studied 
(Fig. I6b[) . It should be noted that the asymmetry of a 
planar Sierpinski gasket fractal could trap the spin in the 
bottom from the spin Hall effect. The up and down spin- 
polarized electrons accumulate on the wide and narrow 
edges of a planar Sierpinski gasket fractal and thus the 
up and down spin electrons can then be selected from the 
asymmetrical shape. Under the appropriate condition, 
there is only one kind of spin electron can pass through 
the asymmetrical device. This geometrical selection pro- 
vides a feasible way of producing pure spin current. 

The transmissions of four different generations are cal- 
culated numerically. As reducing the sites of the devices, 
the widths connection with two leads also decrease. From 
the analysis of different generation of SPG, it was found 
that the contact widths of the leads determine the upper 
bound of the transmission. The transmission affected by 
the connecting width with two leads and the threshold 
of the transmission depends on the number of defects. 
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